import vtk

reader=vtk.vtkMetaImageReader()
reader.SetFileName('C:/Users/bryan/ORTHO/SpineData/Spine2009/CT/C-Spine/C-Spine.mhd')
reader.Update()

marching=vtk.vtkImageMarchingCubes()
marching.SetInput(reader.GetOutput())
#marching.SetValue(0,55)
marching.SetValue(1,105)

mapper1=vtk.vtkPolyDataMapper()
mapper1.SetInput(marching.GetOutput())

actor1=vtk.vtkActor()
actor1.SetMapper(mapper1)


ren=vtk.vtkRenderer()
rw=vtk.vtkRenderWindow()
iren=vtk.vtkRenderWindowInteractor()
rw.AddRenderer(ren)
iren.SetRenderWindow(rw)
ren.AddActor(actor1)

#prop1=actor1.GetProperty()
#prop1.SetColor(.2,.4,.2)


rw.Render()
iren.Start()